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Abstract: Using the AdS/CFT correspondence in the supergravity approximation, we compute 
the energy density radiated by a heavy quark undergoing some arbitrary motion in the vacuum of 
the strongly coupled A/" = 4 supersymmetric Yang-Mills theory. We find that this energy is fully 
generated via backreaction from the near-boundary endpoint of the dual string attached to the 
heavy quark. Because of that, the energy distribution shows the same space-time localization as 
the classical radiation that would be produced by the heavy quark at weak coupling. We believe 
that this and some other unnatural features of our result (like its anisotropy and the presence 
of regions with negative energy density) are artifacts of the supergravity approximation, which 
will be corrected after including string fluctuations. For the case where the quark trajectory 
is bounded, we also compute the radiated power, by integrating the energy density over the 
surface of a sphere at infinity. For sufficiently large times, we find agreement with a previous 
calculation by Mikhailov [hep-th/0305196] . 
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1 Introduction 

One of the basic problems that one can think of in the context of any gauge theory, and in 
particular within a strongly-coupled, conformal, field theory as described by the AdS/CFT 
correspondence [1-3], is that of the radiation by a moving, classical, charged particle. By 
'classical' we mean a particle which is heavy enough to be treated as pointlike and assumed to 
follow a well-identified, classical, trajectory (say, under the action of an external force). And 
by 'radiation' we mean the emission of quanta of the underlying gauge theory which escape at 
infinity, thus generating energy loss. For asymptotically weak coupling, these quanta need to 
be strictly on-shell and hence propagate at the speed of light (for the radiation in the vacuum). 
But in general, the emitted quanta can be also off-shell, in which case they are subjected to 
further evolution (e.g., time-like quanta can decay). In particular, when the coupling is strong, 
we expect such off-shell effects to be very important and generate a very different space-time 
pattern for the radiated energy as compared to weak coupling. 

Consider, for instance, the radiation produced by a heavy quark subjected to a kick, i.e. an 
external force which is localized in space and time. In a classical calculation, which is the same 
as the weak coupling limit of the corresponding field theory, the radiation will propagate away 
from the quark as a spherical shell which is radially expanding at the speed of light (r = t) , with 
a small width Ar determined by the duration of the original perturbation. This is quite different 
from the picture we would expect at strong coupling [4]. There, the radiation should typically 
proceed via the emission of a few virtual quanta, which will then decay into other quanta, thus 
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eventually generating a system of partons with a wide distribution in virtualities. Since time- 
like quanta propagate slower than the speed of light, the energy taken away by those quanta 
should exhibit a rather broad distribution in r at r < t. Since moreover the various quanta 
can be randomly emitted along any direction (in the quark rest frame), this picture also implies 
that the energy distribution should be isotropic (up to a Lorentz boost). This last prediction 
has been checked via an explicit calculation within AdS/CFT of the angular distribution of the 
energy produced by the decay of a time-like wavepacket at strong coupling [5] . 

In view of the above, it appeared as a surprise when other AdS/CFT calculations [6, 7], 
which have also investigated the radial distribution, found that there is no broadening (at least, 
within the limits of the respective calculations): the energy radiated within the vacuum of the 
= 4 supersymmetric Yang-Mills (SYM) theory at strong coupling appears to be as sharply 
localized in r as the corresponding classical result. This was first noticed in Ref. [6] for the 
example of the synchrotron radiation produced by a heavy quark in uniform rotation and then 
extended in Ref. [7] to other situations, including the two problems alluded to above — a heavy 
quark perturbed by a kick and the decay of a time-like 'photon'. As pointed out in [7] (but 
already visible at the level of the calculations in [6]), this lack of radial broadening is to be 
attributed to the fact that, within the supergravity approximation used in these calculations, 
the whole contribution to radiation is generated via backreaction from points near the Minkowski 
boundary of AdSs. 

At this level, it is useful to recall that the supergravity approximation is the classical limit 
of the dual string theory, which neglects string loops and string fluctuations, and is generally 
accepted to faithfully describe the strong 't Hooft coupling limit X = g"^ Nc ^ oo with fixed g <^ I 
of the A/" = 4 SYM theory [1-3]. {g is the Yang-Mills coupling and Nc the number of colors.) 
Furthermore, the 'backreaction' refers to the AdS/CFT calculation of the energy density in the 
gauge theory, which involves the response of the AdSs metric to the 5D stress tensor of the 
bulk object dual to the physical excitation on the boundary. (For instance, this bulk object is a 
Nambu-Goto string in the case of a heavy quark in the fundamental representation of SU(A^c)) 
and a supergravity vector field wave-packet in the case of a time-like photon.) 

A priori, the calculation of the backreaction involves an integral over all the points within 
the support of the bulk stress tensor, say along the string in the case of a heavy quark. Similar 
calculations at finite temperature [8-11] have shown that, in that context, all the points along the 
string provide non-trivial contributions to the energy density on the boundary, which therefore 
shows broadening: points of the string which lie further and further away from the boundary 
provide contributions which are more and more spread in space-time. This is in the spirit of 
the ultraviolet/infrared correspondence [12, 13], which associates increasing distance from the 
boundary of AdSs with increasing virtuality in the original gauge theory. However, in the case 
of the radiation in the vacuum (say, as produced by an accelerated quark), the calculations in 
Refs. [6, 7] show that the integral expressing the backreaction reduces to a boundary contribution 
from the string endpoint at the Minkowski boundary. Thus, there is effectively no virtuality 
involved in this calculation, which 'explains' why the final result shows no spreading. But this 
'explanation' leaves us with a physical paradox, namely why should radiation in a quantum field 
theory at infinitely strong coupling involve only on-shell (light-like) modes, without any trace 
of virtual quantum fluctuations. 
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In Ref. [7] we have also proposed a possible solution to this puzzle, by identifying a class 
of stringy corrections which are not suppressed in the strong coupling limit and which when 
included in the calculation of the backreaction — in an admittedly heuristic way, by lack of 
a proper treatment of string fluctuations in AdSs — seem to provide energy broadening, in 
conformity with the UV/IR correspondence. It would be of course very interesting to make 
further progress with the understanding of string corrections in AdSs, which is an outstanding 
open problem. This is however not the purpose of this paper. Rather, here we shall be more 
modest and extend the results in Ref. [7] in a different direction: we shall provide an exact result 
for the energy density radiated in the supergravity approximation by heavy quark undergoing 
some arbitrary motion in the vacuum of the A/" = 4 SYM theory. 

In spite of our own criticism of the supergravity approximation for the type of problems at 
hand, we believe that the present results are nevertheless interesting for several reasons. First, a 
precise knowledge of the classical, supergravity, result is a first and mandatory step in any effort 
aiming at including string corrections. Second, by itself, this classical calculation is rather non- 
trivial, as it requires an exact, analytic, solution to the problem of the backreaction. Previously, 
such analytic solutions have been given only for particular cases — the most non-trivial one 
being the calculation of the synchrotron radiation in Ref. [6]. Our general results below will allow 
us to simply recover such previous results and extend them to an arbitrary motion. Third, we 
shall explicitly verify that, also in the general case, the whole contribution to the backreaction 
is still coming from the string endpoint near the boundary; hence, in this approximation, the 
radiation propagates at the speed of light, like in a classical field theory. Fourth, our results 
exhibit some other surprising features (besides the lack of radial broadening), which in our 
opinion reflect the limitations of the supergravity approximation: the energy density appears to 
be anisotropic and also negative in some regions of space-time. The anisotropy is unnatural at 
strong coupling, for reasons explained before; it is even more so in the context of the A/" = 4 
SYM theory, where, as we shall see, already the corresponding classical result is fully isotropic^ 
(up to boost effects). As for the negativity of the energy density, which was already noticed (for 
the example of uniform rotation) in [6] , this is in principle acceptable in a quantum field theory, 
but we find it very unnatural in the context of radiation, for reasons to be discussed in Sect. 8. 

Our analysis of the backreaction, that we now outline, will build upon previous constructions 
in the literature. An essential ingredient in that sense is the exact solution, due to Mikhailov 
[15] , for the string profile corresponding to an arbitrary motion of the heavy quark. This solution 
is not fully explicit — it still depends upon a 'retardation time', determined as the solution of 
a generally transcendental equation (see Sect. 2 for details) — , but this is not less explicit than 
the usual textbook treatment of radiation in classical electrodynamics [14], where the results 
are written as a function of the 'retardation time' (the time of emission, related to the time and 
point of measurement by the condition of propagation at the speed of light). By studying the 
energy carried by the accelerated string, Mikhailov has also deduced a formula for the radiated 
power, which appears to be similar to Lienard formula in classical electrodynamics. His results 
have been extended in Refs. [16, 17], where the total energy of the moving quark (proper energy 
plus radiation) has been inferred via a world-sheet analysis. 

^To better appreciate the non-trivial character of this property, one should recall that it does not hold for 
radiation in classical electrodynamics [14]. See also the discussion in Sect. 7. 
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Furthermore, in computing the backreaction, we shall use the general formulas established 
in Ref. [6] which express the energy density on the boundary as a convolution between the bulk 
stress tensor of the string and the bulk-to-boundary propagator. Using Mikhailov's solution 
for the string profile, we shall express the bulk stress tensor in terms of the quark motion on 
the boundary (in Sect. 3), and this will allow us to explicitly perform the integrals yielding 
the backreaction (in Sect. 4). We shall thus find that, due to remarkable cancelations between 
contributions arising from various components of the stress tensor, the only terms which are left 
in the final result for the energy density are boundary contributions from the string endpoint at 
the heavy quark. Then, in Sect. 5, we shall extract the radiative energy density, defined as the 
component of the total energy showing the slowest decay (~ at large distances. This is 

the main result of our paper. By integrating this result over the surface of a sphere at infinity 
(an operation which is well defined when the quark trajectory is confined to a bounded region in 
space), we shall also compute the radiated power (still in Sect. 5). We shall thus find the term 
originally obtained by Mikhailov [15] and also a second term, which is however subleading at 
large times^. In the remaining part of the paper, we shall further discuss our results, compare 
them to some known limits in the literature (in Sect. 6) and also to the corresponding classical 
results (that we shall derive in the context of the A/" = 4 SYM theory in Sect. 7). In our final 
discussion in Sect. 8, we shall emphasize some peculiar features of these results, which point 
towards limitations of the supergravity approximation. 

2 The string profile 

We consider a heavy quark in the fundamental representation of the gauge group SU(A'c) which 
undergoes some arbitrary motion, with trajectory r = rq(t), within the vacuum of the A/" = 4 
SYM theory at strong coupling. (We assume the quark to be arbitrarily heavy, so that the 
notion of classical trajectory makes indeed sense for it.) The quark can either be in isolation, 
or it can be a part of a quark-antiquark pair. The dual supergravity description of the quark 
dynamics is a Nambu-Goto string hanging in AdSs, with one endpoint attached to a D7-brane. 
The string can be either finite, with the other endpoint attached to the same D7-brane (the 
case of a quark-antiquark pair), or it can extend all the way to the center of AdSs (the case of 
a single quark). 

We shall parameterize the AdSs space-time using Poincare coordinates, with metric 
ds^ = Gmn dx^dx"" = ^ {-dt" + dr2 + dz^) . (2.1) 

Here, = (x'^, z) where rr^ = (t, r) are the Minkowski coordinates in the physical space-time 
and z (with < z < oo) is the fifth dimension, also known as the 'radial coordinate in AdSs' 
(not to be confused with the physical radius r = |r|). In these coordinates, the Minkowski 
boundary lies at z = 0, the center of AdSs is at 2; —t- 00, and the D7-brane ends at a distance 
Zm = \/A/ (27rmg) from the boundary, with ruq the quark mass. In what follows we shall choose 
triq to be large enough for Zm to be much smaller than any other interesting space-time scale. 

^Interestingly, this second term is similar to a piece of the total energy of the accelerated string which in 
Refs. [16, 17] has been interpreted as a part of the quark proper (or kinetic) energy. See the discussion in Sect. 6. 
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The string dynamics is encoded in the Nambu~Goto action 
S = -To [ drda^, gab = GuNdaX'^' d^X^ , (2.2) 



where Tq = \f\l2-Kl? is the string tension, r and a are the two coordinates on the world-sheet, 
X'^(r, o") are the string coordinates in AdSs, and gah^ with a, b = t, a, is the induced metric on 
the string world-sheet. 

Choosing t = t and a = z as the two coordinates parametrizing the world-sheet we can 
write the embedding functions and its derivatives as 

xA^ = (t,r„^), X^^ = (l,r„0), xA^' = (0,r;,l), (2.3) 

where a dot or a prime on Vg denote a derivative with respect to t or z respectively. The 
individual components and the determinant of the induced metric read 

g^^ = X-X = |Goo|(-l + r2), <7,, = X'-X' = |Goo|(l + rf), 

gra = X-X' = \Go^\r,-r',, ^ = |Goo| ^/l - f1 + - (r, x r^^, (2.4) 

with |Goo| = L"^ / z"^ . The condition that the action (2.2) be stationary under small variations 
— 7- + 6rs{t, z) yields the string equations of motion 

d_ (l + rf)r,-(r,-rOr^ _ J_ l_ \G o,\[{l - fl)r'^ + {is- r',)rs] ^ ^ 
dt y/1 - r1 + rf - (r, x r',Y |Goo| dz y^l - r2 + _ (r.xr^^ 

As shown by Mikhailov [15], the general solution rs{t,z) to the above equation is implicitly 
determined by 

t = tq+-iqZ, Ts = rq + Vq{t -tq) , (2.6) 

where r^, Vq and 7^ are evaluated at tq, with Vq = Vq the quark velocity and 7^ = 1/(1 — t;^)^/^. 
These two equations should be understood as follows: by solving the first equation (2.6), one 
obtains tq as a function of t and z, which is then inserted into the second equation (2.6) to 
obtain the function rs{t,z). The solutions thus obtained must be restricted to 2: > z^- By 
assumption, Zm is arbitrarily small, but the limit Zm ^ can be taken only after performing 
the 'ultraviolet renormalization', i.e., after absorbing a would-be divergent contribution in that 
limit in the definition of the quark mass. To clarify the physical interpretation of eqs. (2.6), 
notice that they imply 

{rs - Vqf + z"^ = {t - tq)^ and tq{t,z = 0)=t, rs{t,z = 0) =rq{t), (2.7) 

and that the velocity of light in AdSs is equal to one within the present conventions. Hence 
eq. (2.7) can be interpreted as follows: a light signal emitted at time tq at the point with r = Vq 
and 2 = reaches the string (at the point with coordinates Vg and z) at the later time t. Thus, 
eqs. (2.6) show how the string gets built within the bulk via radiation from the quark on the 
boundary. 

We shall often need the derivatives of tq w.r.t. t and z, which read (below, aq = Vq) 
dtq _ 1 dtq _ 7g 



dt 1 + Vq-aq' dz 1 + Z7I Vq ■ ttg 



(2. 
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Using these formulae together with eqs. (2.6) it is straightforward to express the derivatives of 
rg in terms of the boundary motion: 



Vq + 



7g 



zar, 



IqZaq 



1 + Z'J^ Vq-aq 



1 + Vq ■ aq ' 



The following identity is also useful: 



1^ 

7g dt 



(2.9) 



(2.10) 



3 5D bulk stress tensor and backreaction 



The calculation of the space-time distribution of the energy produced by the heavy quark 
requires to solve the 'backreaction problem', that is, to compute the perturbation 5Gmn of the 
AdSs metric associated with the string. For the problem at hand, this perturbation is relatively 
small, of 0{1/N^), and can be computed by solving the linearized Einstein equation with tMN 
(the string stress tensor) as a source. The expectation value (T^^) of the physical stress tensor 
in the boundary gauge theory is then obtained from the near-boundary (z — )■ 0) behaviour of 
5G^i,. Thus, it becomes apparent that we first need to compute the 5D bulk tensor and since 
our source is a string it will be proportional to 6^'^\r — r^). More precisely 

To 



t''^{t,r,z) 



-G 



where it is not hard to see that tMN is given by 

TMN ^0 



r' 



9^^X''x''+grrX''''x'''-gr^ 



(3.1) 



(3.2) 



Now we calculate the components of the above and at the same time lower the indices. The 
metric is diagonal and we do this by multiplying each component by ib|GooP, where we use the 
minus sign only when one of the two indices is equal to 0. After substitution of the common 
coefficient (notice that there is a third |Goo| factor coming from gab, eq. (2.4)) 

TblGooP VX z 1 ^z-fq dt 



-G 27r L3 ^1 - r2 + rf - (r, x r'^f 2-k dtg' 



we have 

^00 

to5 
tij 

ti5 
*55 



(1 + r 



/2^ 

s 



\/A ZJq dt 

27r L3 dtq 

"v/A ZJq dt 
'T^U'dTq 
VX ZJq dt . , 

'2^Udrq'''''''' 

27r L3 dt -'^ ^ 1 + Vs'1^s)[^s^s ~^^s^s)\i 

£7g 5^ 2^ 
27T L3 dtq ^ 



(3.3) 

(3.4) 
(3.5) 
(3.6) 
(3.7) 
(3.8) 
(3.9) 
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Expressing the components in terms of the boundary motion according to eq. (2.9), we finally 
deduce 



^00 



1^ ^ ^ {1 + 2z7,^ v,.a, + zhl{a\ - K x a,f\} , (3.10) 



ioi = ^^^{-^7g4-(l + 2z7,S-«.K- A'[a^-(^,xa,)2]z;^}, (3.11) 



^05 



I? ^ ^ {--7,^ - -^,'[«^ - K X a,f\] , (3.12) 



iij = ^^-JJ ^{^7g(^^g4 + -^^4) + (1 + 2z7g3 vq ■ aq)v\v{ 

+ z2^,6[a2 _ X a,)>^i;^}, (3.13) 

ti5 = ^ ^ ^ {^4 + ^7g^g-ag^^5 + " (^^g x agf]v\} , (3.14) 

= I? ^ § j-:^ + - K x«.)^]| • (3-15) 



4 Energy density in the gauge theory 

The details of this analysis relating the energy density £ = (Too) on the boundary to the string 
stress tensor in the bulk can be found in Ref. [6] from which we shall simply borrow the final 
formulae. Namely, one has 

£{t,r)=£Ait,r)+£Bit,r), (4.1) 

where the two contributions are 

£a = / T- Q{t - t)6"{W) [z{2too - ^55) -{t- t>05 + {x- , (4.2) 

= " 0^ - ^l'(2*oo - 2t55 + tn) - 3(x - xYix - xYUj] . (4.3) 

The argument of tMN is {t, r, z) and the quantity 

W = -{t - tf + {r - rf + z^ (4.4) 

is proportional to the 5D invariant distance between the source point in the bulk and the 
measurement point on the boundary. Eqs. (4.2)-(4.3) are essentially convolutions of Imn with 
the graviton bulk-to-boundary propagator. 

The integration over d^r is trivially done using the (^-function of the string stress tensor given 
in eq. (3.1) and then by using eqs. (3.10)-(3.15) we can express the integrands in eqs. (4.2) and 
(4.3) in terms of the boundary motion. Since rq,Vq and are evaluated at tq, the calculation 
simplifies if we change variable from t to tq . Then for the £a term we have 

■^A = ^ / dtqdz6"{Wq + 2jqEz)[Ao{tq) + zA^{tq) + Z^A2{tq)], (4.5) 

with the definitions 

Wq^-{t- tqf + \r- TqW E ^ {t - tq) - V q ■ {v - T q) = , (4.6) 
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and where the coefficients of the polynomial in z in the square bracket of eq. (4.5) are 

^0 = ?>-iq--iqv\^-il'Vq-aq{t - tg + Vg ■ {r - r^)] + 7g ttg • (r - rg), (4.7) 

Ai = 2^^Vg-ag + ^^[al- {vgxag)'^][t-tg + vg-{r-rg)], (4.8) 

A2 = 0. (4.9) 

Note that, given the ^;-dependencies of the tensor components tjv/TV in eqs. (3.10)-(3.15) and the 
other factors of z in the integrand of eq. (4.2), one would a priori expect the integrand of eq. (4.5) 
to contain a polynomial of second order in z. However, in reality this polynomial is linear since, 
as shown in eq. (4.9) above, the coefficient A2 of the quadratic term is identically zero, due to 
rather non-trivial cancelations. For example, for the terms proportional to z'^[ag — {vg x a^)^], 
there are four different contributions yielding a total coefficient 2^g — 1 — 7q — 7^"^^, which indeed 
vanishes. This has important consequences to which we shall shortly return. 

After replacing t — )■ tg as the integration variable, the argument of the 5-function has 
become linear in z and thus it is easier to perform first the corresponding integration. The 
two derivatives in 5" can be taken w.r.t. Wg and pulled outside the ^-integration. Then the 
5-function sets 

(410) 

Since causality requires H > and z takes only non-negative values, it is clear that the above 
result is non-zero only for Wg < 0. Therefore, 



For the first term it is straightforward to compute the two derivatives. In the second term, 
we first take one derivative, then rewrite the second one as d/dWg = {2E)~^d/dtg, and finally 
integrate by parts to obtain 

^A-5/dV(W,)(^ + A^). (4.12) 

Let tr = tr{t,r) denote the value of tg for which Wg{tg) = 0, that is 

t - tr = \r - rg{tr)\. (4.13) 
Writing 6{yVg) = 6{tg — tr)/2E we finally arrive at 
Ai ^/X 1 d Aq 

One should be cautious to treat tr as a symbolic variable: only after the derivative is performed, 
one can replace tr by its actual dependence on t and r via the solution to eq. (4.13). Also, in the 
previous manipulations we have been a little imprecise about the integration limits in tg and z 
after the change of variables and the associated boundary terms. The most interesting case for 
us here will be a situation where the motion keeps going for ever, meaning —00 < t' < t. Then 
eq. (2.6) implies —00 < tg < t and < 2; < 00, and one can easily check that the integration 
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by parts generates no boundary terms. Indeed, at the upper limit t = tg and the constraint 
eq. (4.13) can be satisfied only when r = rq{t), a situation that we shall not consider. Also, at 
the lower limit tg — )• — oo, eq. (4.13) cannot be satisfied for any finite r. Other situations, where 
the integration domain for t' is finite, need to be considered case by case. 

To summarize, the above integral over tg has support only at tg = t^, where Wg = 0, cf. 
eq. (4.12). Via eq. (4.10), this implies that the integration over z receives contributions from the 
endpoint at z = alone. These special properties are the consequence of the above mentioned 
cancelation of the terms proportional to z'^ in the integrand of eq. (4.5). In turn, they imply 
that the final result (4.14) has the same causal structure as the corresponding classical result, 
that is, as the energy density produced by a source with trajectory rg{tq) in a classical field 
theory. Indeed, the condition = is recognized as the classical retardation condition for the 
propagation of a signal at the speed of light. In particular, tr{t,r) coincides with the classical 
'retarded time' — the time tg at which a light signal must be emitted by the source located at 
rq in order to be received at the point r at some latter time t. Furthermore, the space-time 
pattern of the energy in eq. (4.14) must be the same as in the corresponding classical problem, 
since this is entirely fixed by the trajectory of the source together with the condition that the 
signal propagates at the speed of light. In particular, when focusing on the radiation part we 
expect no quantum broadening to emerge. 

For the calculation of Sb we proceed similarly. Expressing the integrand in terms of the 
boundary motion, and after some tedious but straightforward algebra we arrive at 

"^B = ^ / dtgdz6"'{Wg + 2-fgEz)[B0{tg) + ZB^{tg) + Z^B2{tg) + Z^B^itg) + Z^Bi{tg)], (4.15) 

where the coefficients are given by 

^0 = ^ (r- - Tgf + Jg[{r - Tg) X V g]' , (4.16) 



Bi = --Vg-{r - Tg) + 2^l[{r - rg)^Vg]-{{r - rg)x [^l{vg-ag)vg + ag]] , (4.17) 

^2 = ^ IqV] - 2-ll[{r - rg) X Vg\ ■ {Vg X ag) + jllal - {Vg X agf][{r - Vg) X Vgf, (4.18) 

B3 = Bi = 0. (4.19) 

A priori, the integrand can involve a quartic polynomial in z, but in reality this polynomial is 
just quadratic, since the terms proportional to and z^ have exactly canceled among various 
contributions. Thus the integration has similar properties to the one for since we now have 
three derivatives to take. Once again, the integrations over z and tg are fixed by eq. (4.10) and 
respectively the condition Wg = 0, which together imply tg = tr and z = 0. One finally obtains 

^ _ _B2_ _ ^ 1 ^ _Bl_ _ ^ 1 A /^l A (A 9m 

^ 87r2 73^4 16^2 s 7|H3 167r2 Edtr\E dtr -fgE^ J ■ ^ ' 

The same discussion as for eq. (4.14) applies to potential boundary terms. Eqs. (4.14) and (4.20) 
are our final results for the total energy density produced by the heavy quark. It is interesting 
to notice that Bq and Bi are related as Bi = -^{^gBo), so the last two terms in eq. (4.20) can 
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be combined to yield a somewhat simpler expression for '■ 

^ 87r2 73^4 87r2 E dU \-fgE^ dU E J ' ^ ^ 

In the next section, we shall use these results to extract the energy radiated by the heavy 
quark. But before that, let us perform a first, non-trivial check of these formulae by using them 
to recover the known result for the Coulomb energy of a heavy quark which moves at constant 
velocity (and thus it does not radiate). Assuming uniform linear motion with velocity v along 
the X axis, we have = and then the expressions for the coefficients Aj and Bi simplify 
considerably. After simple manipulations, we deduce 



\/A 72(3 



V 



2\ 



and respectively 

f \/A72 (4 - 2v^)[xl + j^x - vtf] + v^j^x - vtf 

^ 67r2 [xl + ^^{x-vtf]^ ■ ^ ' 

which add together to the expected result [8-10] : 

^ {l + v^)xl + {x-vtf 

12^2 [x\ + -i^{x-vtff ■ ^ ' ' 

Note that all the three terms in eq. (4.20) for 8-q contribute to eq. (4.23), while eq. (4.22) receives 
contributions only from the second, derivative, term in eq. (4.14). 



5 Radiated energy and power 

Prom now on we shall focus on the part of the energy density which is radiated. Following the 
standard definition in the literature, we shall identify the radiation as the part of the energy 
density which falls like l/B? (with R = r — Vqitr)) at large distances from the source. For this 
definition to be meaningful, we shall consider only observation points r which are sufficiently far 
away from the position rq{t) of the quark at the observation time t for the dominant contribution 
of the energy density at r to be falling like 1/ E?. (Indeed, if r is relatively close to rq{t), then 
the retardation condition (4.13) allows for solutions with tr — t and rq{tr) — r, and then 
the energy density at r is dominated by the near-field of the heavy quark, i.e. by its Coulomb 
energy, and not by radiation.) 

To make the power counting with respect to 1/R more transparent, it is useful to notice 
that, for tq = tr obeying eq. (4.13), one has S = R{1 — n-Vq), where we have defined n as 
the unit vector along R. Then, by inspection of the expressions in the previous section, one 
can check that, first, the pieces showing the slowest decay at large distances in Eqs. (4.14) and 
(4.20) are those which behave like as expected, and, second, in order to isolate these 

pieces, it is enough to keep the terms in the coefficients Ai which are proportional to R or 
t — tq and the terms in the coefficients Bi which are proportional to R?. By doing that, one 
eventually finds that radiative contributions oc 1/R? to the energy density come from all the 
terms in Eqs. (4.14) and (4.20) which are proportional to either the square of the acceleration. 
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or to its time derivative (also known as the 'jerk'). Thus, the radiation vanishes in the absence 
of acceleration, as expected. However, unlike what would happen in a classical theory, or in a 
weakly coupled theory at leading order, where the radiation involves only terms proportional to 
the square of the acceleration (see e.g. the discussion in Sec. 7 below), in the present calculation 
at strong coupling we also find contributions proportional to the jerk dq. 

In what follows, we shall exhibit all the radiative contributions to the energy density in 
Eqs. (4.14) and (4.20). In turns out that, in view of the subsequent physical discussion and 
also of the comparison with the respective classical results in Sec. 7, it is meaningful to separate 
between two types of such contributions: (i) those generated by terms in the string stress tensor 
which are by themselves proportional to the square of the acceleration, and (ii) those coming 
from the terms in Eqs. (4.14) and (4.20) which involve derivatives w.r.t. tr- 

(i) Contributions proportional to the acceleration squared, more precisely to the structure 
ttg — (vqXaq)'^, and which are originating from the components (3.10)-(3.15) of Imn, are visible 
in eq. (4.8) for Ai and in eq. (4.18) for B2. They contribute to the energy density via the terms 
without derivatives in Eqs. (4.14) and (4.20), and yield 
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(1) 



87r2 

_v5 7, 
87r2 



(VqXaq)^] 1 + n-Vq 



i?2 



[I -n-Vqf 



(nxvn 



B? (1 — n-Vqj 

These two contributions combine to give the following, relatively simple, expression 



gWu r) - ^ 



%[aq 



^2 _ rq)2(l 



n-Vq)^ 



(5.1) 
(5.2) 

(5.3) 



where it is understood that all quantities related to the motion of the quark (r^, Vq, and a^) are 
evaluated at tq = tr{t, r). Note that the other terms in Ai and B2 do not generate contributions 
of order l/ii^ to the energy density. 

(ii) The remaining terms of order 1/i?^ arise from the 3rd and 4th term of from Bq, 
and from the 2nd term of A priori, that is, within the coefficients Ai and Bi, these terms 
depend only upon the quark velocity Vq and are at most linear in the acceleration a^, but after 
taking the derivatives w.r.t. tr in eqs. (4.14) and (4.20), they generate contributions proportional 
to the square of the acceleration, or to its derivative. Defining ^ = 1 — n-Vq, we find 
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(5.4) 



(5.5) 



where we have neglected derivatives acting on R or t — tq since they generate terms which fall 
faster than Performing the derivative on 1/,^ of the last term in Eq. (5.5) (we do not 

differentiate the unit vector n since this would lead again to terms falling faster than • 
that is, we use dt^£, ~ 



''qji 



'(2) 
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(5.6) 
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Thus, adding the two contributions from eqs. (5.4) and (5.6) we are left only with the last term 
in the last equation, which, after preforming the first derivative w.r.t. tr and returning to the 
original variables, is finally rewritten as 



(5.7) 



If one also performs the remaining derivative w.r.t. tr, one finds that all the ensuing terms are 
proportional to either the square or the derivative of the acceleration, as anticipated. 

The radiated energy density S^ad given by the sum of eqs. (5.3) and (5.7) is the main result 
of this paper. The (relatively high) powers oi 1 — n ■ Vq visible in the denominators of these 
expressions have a kinematical origin: they express the angular collimation of the radiation due 
to the Lorentz boost, which was to be expected, independently of the value of the coupling. In 
the ultrarelativistic limit ~ 1 or 7^ 1, one can write (with a denoting the angle between 
the vectors Vq and n): 1 — Vqcosa ~ (l/2)(l/72 + a^). This makes it clear that the radiation 
is emitted within a small angle a ~ l/7g around the direction of the quark velocity Vq{tr), so 
like for the corresponding classical problem [14]. 

Using the above results for £^rad = "^^rad + "^rad' shall now compute the radiated power. 
By integrating the energy conservation law dt{T^^) + di{T^^) = over the whole space and 
using Gauss' theorem together with (T*^*) ~ n^^T^^) for the dominant respective contributions, 
proportional to at large distances, one finds (recall the notation {T^^) = £) 

d F r 

- — = lim / dn8{t,r), (5.8) 

dt r-i>oo J 

where the left hand side represents the energy radiated per unit of observation time. In practice, 
it is more convenient to define the power Prad as the energy radiated per unit of emission time 
tr- Then by using the above formula together with dt/dt^ = sees that the power radiated 

per unit solid angle reads (below, R — )• 00) 

dfiad ^ 
dO dtr 

Note an important, implicit, assumption in the above argument: we have made the hypothesis 
that, at all the points on a sphere at infinity (r — t- 00), we have R=\r — Tg(tr)| — (which in 
turn implies that only the far-zone contributions oc l/B? to the energy density and fiux have to 
be retained). This is correct provided the motion of the quark is hounded, such that its trajectory 
rq{tr) does not cross the sphere at infinity. More precisely, it is enough that this condition be 
satisfied during the acceleration phase of its motion, since this is the only phase which creates 
radiation. Thus, our subsequent results for Pj-ad a-re only valid provided there exists some fixed, 
but arbitrary, distance such that Tq[tq) < tq for any tq within the acceleration phase. We 
shall later make some comments on the case of unbounded motion. 

Some standard and useful integrals to perform the angular integrations in eq. (5.9) are listed 

"rad 



R^Srad Pr.d= dn{l-n-Vq)R''£,,d. (5.9) 



in Appendix A. For the first contribution coming from S^ll we find 



^rS = ^7,'K-Kx«.)']. (5-10) 

which is the result inferred in [15] from a world-sheet analysis (i.e., without an explicit cal- 
culation of the backreaction, but merely via a calculation of the energy flux down the string). 
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Remarkably, this expression has the same structure as the respective classical result, that is, the 
Lienard formula in classical electrodynamics [14], that we shall extend to the case of the A/" = 4 
SYM theory in Sec. 7. For the second contribution from £^J^ we obtain 

The fact that this term is a total derivative w.r.t. the emission time tr rises some puzzles for its 
interpretation as a contribution to the radiated energy (see the discussion at the end of Sect. 6). 
It is therefore interesting to notice that a term with a similar structure has been interpreted in 
Ref. [16, 17] as a contribution to the proper energy of the quark (and not to its radiation). We 
shall return to this point in the next section. 

Eqs. (5.10) and (5.11) are our final results for the radiated power. It is important to keep in 
mind that these results have been derived here for the case of a bounded quark motion. Because 
of that, in evaluating these formulae one can use the approximation tr — t — r for most (but not 
all) purposes. 

6 Applications 

Here we shall apply the general results derived in the previous section for the radiated energy 
density and the power to specific quark motions. 

(i) Uniform rotation: We shall first recover the results for uniform circular motion originally 
obtained in [6]. Using spherical coordinates r = (r, 9, (p) and parametrizing the boundary motion 
as 

rg{tr) = {Ro,TT/2,Ujtr), (6.1) 

our expressions in eq. (5.3) and eq. (5.7) lead to the following two contributions to the density 
of the radiated energy 

(1) _ ^ 

(2) ^ ^00^ A-lj-Av^ smH + ^e 

rad 247r2r2 72^6 ' vo-^j 

with V = loRq, a = lo'^Rq, and where, according to our earlier definition above eq. (5.4), we have 

^ = 1 — V sin 9 sm{(j) — utr) ■ (6-4) 

Note that the contribution in eq. (6.3) is fully arising from the last term, proportional to n aq, 
in eq. (5.7), since Vq-ag = for the problem at hand. 
Adding the two pieces in eqs. (6.2) and (6.3) we find 

_ ^/Xio^ 4 - 7^ - 4t;2 sin^ 9 + 
■^rad - 24^2^2 ^2^6 ' V6.5j 

which is indeed the same as the result for the radiated energy density in [6] (cf. eq. (3.72) there). 
The energy density (6.5) is proportional to 1/^^, hence in the ultrarelativistic limit f ~ 1 or 
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7 S> 1, it is strongly peaked at the minima of ^, defined by sin 6 sm{(f> — cotr) = 1- This condition 
describes a spiral in the plane 6 = vr/2, located at (recall that tr — t — r) 

^{t,r) ~ ^+u{t-r). (6.6) 

Using 1 — vsm.9 ~ (l/2)(l/7^ +cP') where a = tt/2 — 0, it is obvious that the energy is localized 
within a small angle a ~ I/7 around a = 0, or0 = 7r/2. This is the expected collimation due 
to the Lorentz boost, as already discussed in relation with the general formulae (5.3) and (5.7). 
Furthermore, using (6.5) one can check that the arms of the spiral have a tiny radial width^ 
Ar ~ Ro/l^ [6], exactly like in the corresponding classical problem [14]. As explained in the 
Introduction, this feature is very surprising in a quantum theory at strong coupling, where one 
would rather expect broadening due to the virtual quantum fluctuations. In the context of our 
calculation, this follows from the fact that, as explained on Sec. 4, the whole backreaction arises 
from the string endpoint at 2: = 0. (See also [18] for a different perspective of this problem.) 

We now compute the radiated power in terms of the quark's own time. From eq. (5.11), it 
is clear that P^.^^^ = in this case, so the only contribution comes from eq. (5.10) and reads 

Prad = ^7V. (6.7) 

This coincides with the respective result in [6] and also with the result of the world-sheet analysis 
in [15]. 

(ii) Non— uniform circular motion: It is of course straightforward to apply our general 
expressions to an arbitrary circular motion, but to be more precise we shall focus on the specific 
motion 



rg{tr) = iRo,7r/2,^g{tr)) with ^g{tr) = VW+^/Ro, (6.8) 

for which we shall directly compute the radiated power. The (angular) velocity, whose magnitude 
approaches the speed of light at large times, and the acceleration are given by (below and e^p 
are the respective unit vectors) 

-.(t.) = ^^7== and = - ^ + j^^^ e„ (6.9) 

so in particular 7^ = y^t^ + b'^/b and Vg-aq is not vanishing anymore (in contrast to the case of 
uniform rotation). Hence both terms contributing to the radiated power, (5.10) and (5.11), are 
now non-zero, and this is interesting as it allows us to observe an hierarchy among these terms, 
that we believe to be generic. Namely, Pj^^^ dominates over pI:^\ in the ultrarelativistic limit 
7g ^ 1, and hence also for sufficiently large times (in the problems where the velocity grows 
with time, due to acceleration). Specifically, for the motion in eq. (6.8), one finds 

(1) _^ 4 + b^Rl p(2)__^{AJ_ 



■'in order to deduce this property from eq. (6.5), it is not enough to use the simpHfied version of the retardation 
time tr — t — r; rather, one needs a more precise analysis of the retardation condition, which also takes into account 
the colhmation of the radiation along the direction of emission; see [6, 14]. 
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As anticipated, the first term P^^^ is the dominant one for large enough times such as t'^ S> 6i?o- 
This example also gives us some insight into the physical origin of this hierarchy: eq. (5.10) for 
-^rad involves the component of the acceleration which is transverse to the velocity^ (the radial 

(2) 

component of in eq. (6.9)), whereas P^:^^ in eq. (5.11) rather involves the respective longitu- 
dinal component (tangential in the case of eq. (6.9)). We thus recover a feature familiar in the 
context of classical radiation [14]: rotation is much more effective than tangential acceleration in 
producing radiation, since the velocity Vq changes rapidly in direction while the particle rotates, 
even though its change in magnitude is relatively small (or even zero for uniform rotation). 

(iii) Uniform linear acceleration: A classical particle subjected to a constant force F = Fei 
follows a trajectory Xq{tr) = yt^ + fe^ where x = x^, b = m/F and we selected convenient initial 
conditions at t = 0. Clearly, this motion is unbounded: the trajectory will eventually cross the 
sphere 'at infinity' that we use to define the radiated power. In view of that, we do not expect 
our previous results for the power to also cover this case. Notwithstanding, let us first see what 
these results yield if naively applied to this case. Using 



and hence ^^VqQq = trjl? ^ one easily finds that the two contributions in eq. (5.10) and respec- 
tively eq. (5.11) are now of the same parametric order and thus contribute on the same footing 
to the final result for the power, in contradiction with our general expectations and also with 
the previous examples. Namely, one (naively) has 



(1) _ ^ (2) _ ^ p ^ .2) 



Moreover, this result also contradicts the independent calculation in Refs. [10, 19], where the 
radiated power has been extracted from a world-sheet analysis, as the energy flow across the 
induced horizon at 2; = 6. That previous calculation furnished a result equal to P'^^ in the above 
equation, which if course would be also the prediction of Mikhailov's analysis for the problem at 
hand [15] (since in that analysis the total power reduces to our P'^^- Clearly, these mismatches 
shed further doubts on the validity of the above calculation, that is anyway transgressing the 
validity limits of our general calculation. 

Let us therefore redo our analysis of this problem, but in such a way to stay within the limits 
of our general discussion. To that aim, we assume that the quark is under uniform acceleration 
only for a finite period of time toi £^nd we measure the radiated energy at a distance r S> to- 
By doing this we effectively reduce the motion to a bounded one. Let us be more specific and 
consider the one-dimensional motion 



t^tr + I) 



Xqitr) = e{-tr) b + e(t,)G(to " tr)VW+^+ ©(^r " ^o) (6.13) 

where can be taken to be much larger than b so that the quark becomes eventually ultra- 
relativistic. (The last term in eq. (6.13) describes a constant velocity motion with the velocity 



More precisely, the term in eq. (5.10) receives contributions from both the radial and the azimuthal 
components of the acceleration in eq. (6.9), but the dominant contribution at large times, represented by the term 
proportional to tf. in the numerator of Pj^'j in eq. (6.10), is generated by the radial piece of a,. 
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acquired at tr = to-) In this last example we shall evaluate the total energy radiated, and 
therefore we need to integrated the total power over tr- Then it is straightforward to see that 
eq. (5.11) will not contribute to the final results, since it involves a total derivative w.r.t. and 
the acceleration vanishes outside the interval [0, to]- Thus the power in eq. (5.10) will determine 
the total energy radiated which is 

E„.^§'l,. (0.14) 

in agreement with [10, 19]. 

(2) 

This last calculation also illustrates a rather curious feature of Pj^.^^ in eq. (5.11), which 
makes us feel uncomfortable about its physical interpretation as radiation: this term yields no 
contribution to the radiated energy for any motion where the acceleration is non-zero over only 
a finite interval of time, and also for any periodic motion for which one computes the total 
radiation over one period, or an integer multiple of it. As already mentioned after eq. (5.11), 
a term with a similar structure appears in the world-sheet calculation of the total energy of 
the moving quark (which is the same as the total energy carried by the string) in Ref. [16, 17]. 
Specifically, by taking the limit of a very heavy quark [zm — ^ 0) in eqs. (2.32)-(2.33) of Ref. [17], 
one obtains the following expression for the quark energy 

Eq{t) = mg-fq{t) - ^j^Vg-ag + y dig 7^ [a^ - (-u^ X a<;)2] , (6.15) 

where the first two terms in the r.h.s. are interpreted [16, 17] as the quark proper (or kinetic) 
energy, while the third one, which is the time-integral of P^l^ in eq. (5.10), as the radiation. As 
anticipated, the second term in the above equation has the same structure as the contribution to 

(2) 

the 'radiated energy' that would be obtained from P,^.^^, eq. (5.11), but with a different numerical 
coefficient. Recall that the reason why we have identified this term as radiation in Sect. 5 was 
because it arises from a piece in the energy density which falls off like l/i?^ at large R. This 
suggests the interesting possibility that a piece of the quark proper energy have a tail at large 
R which cannot be distinguished from radiation. We leave this question, as well the calculation 
of the total energy via the backreaction, for further studies. 



7 The classical result 

In the previous discussion, we have already anticipated some similarities between the predictions 
of the supergravity approximation for the strong coupling limit and the corresponding results in 
the classical approximation, which are also the leading order results at weak coupling. For this 
comparison to be more precise and in preparation of the physical discussion in the next section, 
in this section we shall explicitly solve the corresponding classical problem — the radiation by 
a heavy quark undergoing some arbitrary motion in the A/" = 4 SYM theory at weak coupling. 
To our knowledge, the general result that we shall derive here has not been presented elsewhere, 
except for the case of the uniform circular motion that was discussed in [6]. But even in that 
case, the final results and the associated physical discussion have been plagued by some mistakes 
in the numerical factors that we shall here correct. 
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The general structure of the classical theory describing a massive test quark^ propagating 
through the vacuum of the A/^ = 4 SYM theory has been clarified in Refs. [6, 20]. As explained 
there, the heavy quark radiates both vector (gauge) fields and scalar fields, and in the limits of 
interest here (arbitrarily weak coupling and very large quark mass) this radiation is described by 
decoupled, linear equations, which generalize Maxwell equations to the theory at hand. These 
equations are then solved in the standard way, to give 



47r(l - n-Vg)R 



{l,Vq) and X 



47r7g(l -n-Vg)R ' 



(7.1) 



with A^^ the vector field, x t^is scalar field, R = r — Vg and n the unit vector along R as earlier, 
and the proper counting of the color degrees of freedom for the radiated field is encoded in 
Cgjj = A/2. As before, the above expressions are to be evaluated at the retarded time ty which 
is the solution to eq. (4.13). In general the energy density is obtained from 
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with B the magnetic field. Since we are interested in the radiated energy, we keep only the 
contributions which fall like 1/R when computing the electric field and the derivative with 
respect to time of the scalar field. This yields 
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|(Vx)rad| for the radiation, we deduce that 
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(7.5) 
(7.6) 



where we have substituted e^g 
n ag cancel^ and we obtain a simple result: 



A/2. Adding the two contributions the terms depending on 



cclass 
^rad 



A - Kxa,)^] 



327r2 (^r-rg)'^{l-n-Vg)^ 



(7.7) 



It is very interesting to notice that with the replacement A — >• 4\/A, this is exactly the same as 
the £^l\ piece of the strong coupling result in eq. (5.3). Since by assumption r ^ Vg and thus 
the retarded time can be approximated as t,. ~ t — r, we see that the only angular dependence is 



2 hypermultiplet, that is in the funda- 



^More precisely an infinitely massive spin-1/2 particle from the A/" 
mental representation of the SU(A'^c) gauge group. 

^This cancelation for the particular case of uniform circular motion was not realized by the authors of [6], 
because their corresponding expressions (2.20a) and (2.20b) for ^vector and ^scalar miss a numerical factor of 1/2 
and 2 respectively. 
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the boost factor (1 — n-Vq)^ in the denominator. This means that in the non-relativistic hmit 
eq. (7.7) becomes isotropic and reduces to 

1 A 

^'^'^ - 327r2 r2- ^'■^> 

Obviously this is a property which is not shared by QED, where only vector fields are radiated, 
and the radiated energy as given in eq. (7.5) contains anisotropic pieces. 

One can compute separately the vector and scalar contributions to the power (see again 
Appendix A for the corresponding integrals), which read 



-Pvector = 7g [ttg - (^g X flg)^] and Pgcalar = ^ 7g ["g " (l^g X a,)^] , (7.9) 

leading to a total power 



^r'^ir = ^ 7,' [a^ - (^, X a,)2] , (7.10) 

which, up to the replacement A — )• 4\/A, is the same as the piece pj^^] of the corresponding 
supergravity result, cf. eq. (5.10). 



8 Discussion and open issues 

One of the main results of this paper is the verification of a conjecture made in [7] that for an 
arbitrary relativistic motion of a heavy quark, and in the supergravity approximation to the dual 
string theory, it is only the endpoint of the dual string at z = which contributes to the radiated 
energy. That was first observed in [6] for the uniform circular motion and then extended in [7] 
to a general non-relativistic motion and also, mutatis mutandis, to other types of radiation, like 
the decay of a time-like wave-packet. (The dual description of the time-like wave-packet is 
a supergravity field falling into AdSs. Then the corresponding statement is that the radiation 
is generated only from the starting point of the trajectory at z = 0.) This property implies 
that the radiation propagates at the speed of light and therefore the space-time distribution 
of this radiated energy is very similar to that of the corresponding classical radiation, without 
any sign of quantum broadening [6, 7]. As argued in [7] this leads to a radial distribution 
which is to difficult to reconcile with quantum mechanics, which sheds doubts on the validity 
of the supergravity approximation as the correct, dual, description of the strong-coupling limit. 
Moreover, it was shown there, for a specific example and via an admittedly heuristic calculation, 
that there are particular string corrections (which in the light-cone gauge appear as fluctuations 
in the longitudinal coordinates of the string points) which are not suppressed when A — >■ oo, 
and hence should be treated as a part of the leading order theory at strong coupling. One effect 
of those fluctuations (at least within the limits of the calculation in [7]) is to provide a radial 
broadening for the energy distribution, in agreement with expectations from both quantum 
mechanics and the UV/IR correspondence. 

The proper calculation of string fluctuations in a curved space-time is an outstanding open 
problem, that we shall not attempt to address here. Rather, we would like to emphasize some 
curious features of the previous results obtained in the supergravity approximation, which look 
rather implausible to us on physical grounds and may represent additional shortcomings of this 
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approximation (besides the lack of quantum broadening). The peculiarities to be discussed here 
are all associated with the second contribution to the energy density, S^:^!^ in eq. (5.7). 



(i) The lack of isotropy in the non-relativistic limit 

As mentioned in the Introduction, at strong coupling one expects the radiation to be isotrop- 
ically distributed at large distances away from the source [4, 5], except for the trivial anisotropy 
introduced by the Lorentz boost. Indeed, the radiation should typically proceed via the emission 
of off-shell quanta'', which then should evacuate their virtuality via successive branchings. This 
gives rise to a partonic cascade through which the original energy and momentum get divided 
among many quanta. Due to their large number and to the absence of any preferred pattern in 
the process of branching (at strong coupling), these quanta should have an isotropic distribution. 
This is generally not the case at weak coupling (say, in classical electrodynamics), although it 
happens to be the case in the weak coupling limit of the = 4 SYM theory, as shown in Sect. 7 
(see also [5]), because of the additional symmetries of this theory. 

In view of the above, we find it extremely surprising that the supergravity result for the 
radiated energy density is not isotropic in the non-relativistic limit, especially in the context of 
the M = 4 SYM theory, where the isotropy is realized already at weak coupling. Indeed, when 
Vq <^ 1, Eqs. (5.3) and (5.7) reduce to 



which is manifestly not isotropic. Moreover, this last, anisotropic, term can even dominate 
over the first one in some cases, as shown by the example of the uniform rotation: there, 
a = lo'^Rq = Lov and a = io'^ Rq = uj'^v, so clearly d ^ when f <C 1. Since moreover the sign of 
n dq is oscillating when changing the direction of observation, one sees that the radiated energy 
density is negative in some regions, which brings us to our second puzzle. 

(ii) The negativity of the radiated energy density 

In Ref. [6] already the authors noticed that the energy density (6.5) radiated in the case 
of uniform rotation can become negative in some regions of space-time. From our present 
discussion we know that this behaviour must be associated with the second piece (6.3) of the 
radiation, which in some regions can become negative and also larger in magnitude than the first 
piece (6.2). In Ref. [6], where only the relativistic case was considered, the regions of negative 
energy were relatively small (and localized near the edges of the arm of the spiral) and besides 
the negative values reached by the energy in those regions were much smaller than its positive 
values towards the middle of the spiral arm. In the non-relativistic case, however, we have just 
seen that (for the case of rotation at least), the second piece (8.2) dominates over the first one 

^This follows from the uncertainty principle: quanta with a large virtuality Q have a short formation time 
At ~ LojQ^ ^ where lo is the energy carried by the quanta. However, after being emitted, such quanta need to 
further radiate to become on-shell, which explains why their emission is suppressed at weak coupling. 



,(1) ^ a| 

rad g^2 ^2 




which is isotropic and similar in structure to eq. (7.8), and respectively 
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(8.1) anywhere except at the particular points where n dg vanishes, and that the sign of this 
second piece oscillates. Specifically, the non-relativistic limit of (6.5) reads 



which arises, as expected, from the n dg term of £^,^1^ in eq. (8.2). Clearly, this (dominant) 
contribution to the radiated energy density oscillates around zero with the positive and negative 
maxima being of equal magnitude. 

In principle, regions of negative energy density can occur in a quantum field theory, in the 
process of subtraction of the ultraviolet (UV) divergences. For the problem under consideration, 
such UV issues could affect the proper energy of the heavy quark, as carried by its near field, but 
on the other hand we find them rather unnatural in relation with the far fields and the radiation. 
As manifest say in eq. (8.3), the length scale associated with such space-time variations in the 
radiated energy is not some UV cutoff, but rather is determined by the external force that is 
giving the quark the specified motion. 

Note also that eq. (8.3), and more generally the n dg term of eq. (8.2), do not contribute 
to the radiated power, since they integrate to zero. The power appears to be dominated by the 
first term (5.10), and thus be positive, for all the examples that we have investigated. 

To summarize, the anisotropy and the negativity of the energy density associated with 

(2) 

the contribution £^^^ in eq. (5.7) look very unnatural to us and make us feel skeptical about 
this particular term. In our opinion, these unappealing features are merely an artifact of the 
supergravity approximation which will be corrected after including string fiuctuations. It is also 
possible that this term, or at least a part of it, represent the tail of the quark proper energy at 
large distances, as suggested by the comparison between the associated power, eq. (5.11), and 
the results in [16, 17] (cf. the discussion after eq. (6.15)). 

Also, the fact that the problems alluded to above are solely generated by the second piece, 
eq. (5.7), of the energy density does not mean that we fully trust the other piece in eq. (5.3). 
In spite of its rather appealing structure and of its similarity with the corresponding classical 
result, this term too has been produced from the string endpoint at 2; = and thus it shows no 
quantum (radial) broadening. We therefore believe that also this term will be modified by string 
fiuctuations, in the sense of acquiring a spread, but in such a way that its spatial integral giving 
the power will remain unchanged. Indeed, we believe that the correct result for the power (at 
least for a bounded motion and sufficiently large times) is given by P^f^^] in eq. (5.10), because this 
expression has been suggested by independent considerations (based on a world-sheet analysis) 
in Refs. [15-17] and because it coincides with the energy flow at the world-sheet horizon in all 
the examples that have been worked out in the literature. 

At this point, we should recall that Ref. [5] has studied the angular distribution of the 
energy density produced by the decay of a time-like wave-packet within AdS/CFT, and found 
that this is isotropic in the supergravity approximation and it is only weakly affected by string 
fluctuations (at least, in a heuristic treatment of the latter inspired by flat -space string quanti- 
zation). However, in that particular problem there was no kinematical scale which could induce 
an anisotropy (the wave-packet was spherically symmetric and at rest), unlike in the heavy 
quark problem under present consideration, where there are such scales. For us, our present 
results signal that, in general, the supergravity predictions cannot be trusted neither for the 




(8.3) 
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angular distribution of the radiation, nor for the radial one. Therefore, any progress towards 
better understanding the effects of the string fluctuations would be of paramount importance. 
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A Useful integrals 

Here we list some standard integrals which are useful when calculating the total power. With v 
the velocity, 7 the Lorentz boost factor, a an arbitrary vector and n the unit vector 

n = (sin^cos(/), sin ^ sine/), cos (A.l) 

we have 



/ 



dn—^— = 4n^\ (A.2) 
(1 — n-vy 

,^ na 167r r . 
and also 

dn = 87r72 tanh-i v, A.5 

(1 — n-vY V 

dn Z"'''''^' = 47r7' - — tanh"! v. (A.6) 
(1 — n-vy V 

An easy way to perform all the above or similar integrations is to assume, without any loss of 
generality, that instantaneously the particle is moving along the third axis, that \s v = (0, 0,u). 
Then it is straightforward to perform the integral / dil [[n-aY / [a — v cos 6)] with integer p >0 
and arbitrary a > v. The desired integrals follow by an appropriate number of differentiations 
with respect to a evaluated at a = 1. 
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